Gravitational lens equation: critical solutions 
and magnification near folds and cusps 
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We study approximate solutions of the gravitational lens equation and corresponding lens magnification factor 
near the critical point. This consideration is based on the Taylor expansion of the lens potential in powers of 
coordinates and an introduction of a proximity parameter characterising the closeness of a point source to the 
caustic. Second-order corrections to known approximate solutions and magnification are found in case of a general 
fold point. The first-order corrections near a general cusp are found as well. 
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INTRODUCTION 

Main equation of gravitational lens theory (eq. 
([1]) below) sets a relation between the angular posi- 
tion y of the po int source and the observable position 
X of its image |12| . The main interest is related to 
critical points of the two-dimensional lens mapping, 
i. e. the values of Xcr where Jacobian of the lens map- 
ping vanishes: J (xcr) = D {yi,y2)/D {xi,X2)\^^^ = 
0. In its turn, the image of set of critical points is a 
set of caustic curves. Each caustic typically appears 
as a closed smooth curve with the so-called cusps at 
some isolated points. Regular points of a caustic as 
singularities of differential mapping are called folds. 
When a point source crosses the fold caustic, the two 
critical images either emerge or disappear. The crit- 
ical images of the point source approach the critical 
curve and their brightness tends to infinity when the 
source comes close to fold. In the vicinity of a cusp 
we have three critical images with infinite bright- 
ness, but only two images disappear afte r crossing 
the cusp. The details can be found in [12]. 

The singular properties of caustic points play a 
key role in explanation of some qualitative features 
and quantitative characteristics of the gravitational 
lensing phenomenon. For example, the qualitative 
picture of quadruple lensing can be modelled by a 
singular isothermal ellipsoid [4j. Specifically, rela- 
tive positions of four images and their brightness de- 
pend on position of the source with respect to the 
caustic [21 [T^. Another example is related to the 
so-called strong microlensing events, which are inter- 
preted as a crossing of a microcaustic by an extended 
source. In this case astronomical observations give 
us a chance to get some information about size of 
the source and distribution of the brightness on its 



surface [SEKII]. 

The well known approximate solutions of the lens 
equation and expressions for magnification of each 
image obtained in the lowest approximation [6l 13 [12] 
have a sense of asymptotic relations, which are per- 
formed the better, the closer the source is located 
relative to the caustic. 

In the case of a fold caustic, the first-order cor- 
rections for approximate coordinates of the critical 
images were found in [2j; corrections for the mag- 
nification of separate images were obtained in pa- 
per |1Q) . Note that it is possible to observe only a 
total brightness of all microimages during a strong 
microlensing event. In this case, the first-order cor- 
rections for magnification of two critical images are 
mutually cancelled. The second-order corrections for 
image coordinates, as well as for the magnification, 
were found in papers [UlS]. Besides, it was demon- 
strated on example of the strong microlensing event 
in image C of gravitational lens system Q2237+0305 
that the second-order corrections can be statistically 
significant. It was made under the simplifying as- 
sumption that there is no continuous matter near 
the line of sight [H [3] . Because of importance of ac- 
counting dark matter, we generalize expressions for 
the second- order corrections near the fold in the 
present paper. 

Concerning the cusp caustic, the first-order cor- 
rections were considered in [4J, but some expressions 
in that paper require revisions. Moreover, calcula- 
tions are missing logical conclusions, they were left 
on some intermediate stage. Therefore, the second 
part of our paper is dedicated to looking for com- 
plete and more compendious expressions in the first- 
order approximation for the coordinates of images 
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and magnification near the cusp. 

LENS EQUATIONS NEAR CRITICAL POINT 
The normahzed lens equation has the form: 



y = X - V$ (x) , 



where $ (x) is the lens potential. This equation re- 
lates every point x = {xi,X2) of the image plane 
to the point y = (1/1,2/2) of the source plane. In 
the general case, there are several solutions X(^) (y) 
of the lens equation ([1]) that represent images of a 
point source at y; we denote the solution number by 
the index in parentheses. 

Potential $ (x) satisfies equation A<I> = 2k, where 
k (x) is the density of continuous matter on the line 
of sight normalized on the so-called critical den- 
sity. The magnification factor of each separate im- 



J (X(0 (y)) ' where J (x) 



age is K(^i) (y) = 

\D{y)/D(x)\ is the Jacobian of the lens mapping 
©• 

Recall that, critical curves of mapping ([1]) are de- 
termined with equation J (x) = 0. Caustics are im- 
ages of critical curves obtained with mapping ([1]). 
The stable critical points of a two-dimensional map- 
ping can be folds and cusps only. 

Using standard approach to examine neighbour- 
hood of the caustic, potential near the point p^r 

of 

the critical curve can be approximated with the Tay- 
lor polynomial. Let this point be the coordinate ori- 
gin. We suppose that eq. ([TJ maps Pcr onto the co- 
ordinate origin of the source plane. Then, we rotate 
synchronously the coordinate systems until the ab- 
scissa axis on the source plane becomes tangent to 
the caustic at the origin; the quantity \y2\ defines 
locally the distance to the caustic and yi is a dis- 
placement along the tangent. 

With a sufficient accuracy, the lens equations 
have the following form: 

yi = 2 (1 - ko) xi + aixf - 02X3 + 262a^ia:;2+ 

+ C2X1 — 3ciXiX2 — diX2 + 3d2XiX2 + §1X2 + •■• 



y2 



b2x\ - hix\ - 2a2a;ia;2 -|- d2x\- 



— "idlXlX^ + C2x\ — 'iClX2x\ + /3X2 + ... (2) 

Here ko = k (0) is the matter density at the origin 
and the following notations are: 

ai = -^>,iii/2;a2 = ^>, 122/2; 61 = ^>,222/2; 

62 = -^>, 112/2; Ci = 1122/6; C2 = -^>,2222/6; 
di = $,1222/6; d2 = -$,1112/6; 5 = -$,12222/24; 

/ = -$,22222/24. 



When density k is constant, then oi = 02 = a, 
hi =h2 = b, ci = C2 = c, di = d2 = d. Parameter d2 
will not appear in the following formulae; therefore 
we put di = d. 

APPROXIMATE FORMULAE 

NEAR FOLD CAUSTIC 

One of the approaches for finding critical solu- 
tions of eq. ([1]) involves an expansion of the image 
coordinates into series in powers of some parameter 
t, which demonstrates proximity to the caustic p]- 
3 . If we put yi = t^yi, then, as it was shown in 
l-[3], the critical solutions of eq. ([1]) are analytical 
:"unctions of parameter t, and xi = t^xi, X2 = tx2, 
where xi (t), X2 (t) are zero-order functions. Putting 
these expressions into Taylor expansion of eq. ([TJ, 
and restricting our solutions to second-order terms 
inclusive, we get the following equations: 



2/1 = 2(1- ko) xi - 02^2 + t (262X1X2 - dx2) + 



-|- t"^ (aiXi — 3ciXiX2 -|- 5X2) 



?/2 



-61X2 + t (-202X1X2 -I- C2X2) + 

+ (62X? - Uxixl + fxfj . (3) 



After performing calculations, it is enough to put 
t = 1 and thus return to the initial variables y^ and 

Xi . 

A condition that initial point fold is that 

hi 7^ 0. Without losing generality of our approach, 
we assume that hi < 0. When density k is constant, 
then oi = 02 = a, 61 = 62 = b, ci = C2 = c. There- 
fore, the system Q includes four additional parame- 
ters in comparison with previous case of [T]-[3] where 

k (x) EE 0. 

We seek solutions of equations ([3]) accurate within 
second-order terms in a form: xi = xio + xut + 
Xi2t^, X2 = X20 + X2it + X22i^- Imposing notations 
R'^ = 02 + bib2, (7 = 1 — ko and e = ±1 we find 
the following expressions in the zero-order approxi- 
mation: 



1 



= Trim- 022/2/61) , X20 = e^/y2/\bl\. (4) 

Two signs of parameter e correspond to two critical 
solutions. The first approximation gives: 



Xll 



X21 



[a2R^ - {hid + a2C2)cF\y2] , (5) 
-02612/1 + (al - C2(t) y2 



26fa 



(6) 
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Concerning the second-order approximation for the 
first coordinate we found: 



Mivl + M2yiy2 - M^y^ 



r,2 



Xl2 



(7) 



where 



Ml = h\ (3a26i&2 + 2al - aihl) , (8) 



M2 = 2bi [bl (aias - 2hl - SciO") - 61 (7a|62- 

-(62C2 + 6a2(i)cr) - - C2(t)] , (9) 

M3 = 6f [aia2 - 40262 + (462^ - 60201)0" + Aga"^] + 

+ hi [-110^62 + (16o|d + 6026202)0- - (402/ + 

+6c2d) a^] - 602 [al - C2a)^ . (10) 

And for the second coordinate: 

r——Niy2 + iV2yi + N^yl/m ,^ , , 
^MM ^2 ' (11) 



a;22 = e 



= -5o2i?^ + 10 (o26id + O2C2) a- 

- (5c2 + 4/6i)a2, (12) 

N2 = 661 [a2R^ - {bid + 02C2) a] , (13) 

N3 = -bjR^. (14) 

In its turn, for the Jacobian of the lens mapping, 
calculated in points where images are situated, we 
found: 



J = tjQ + t^Ji + rj2 



Jo = 4:ea^/\h\§2, Ji = 4 



- C2(7 



(15) 

m, (16) 



J2 = e^/y2/\hl\ '- — , (17) 



p = 2^62/61 -r/86f 



3^2 



(21) 



T = 61 [190^62 - 4010261 - (30o2d + I262C2 - 

-I261C1) a + 12/(7^] + 15 {al - C2af , (22) 



Q 



3o2 — 2oi6f + 5026162 — 3 (02C2 + 6id) a 



462fj2 



i?2 



2|6i|(t2 



(23) 
(24) 



In comparison with the formulae that were found un- 
der assumption of /c = 0, we shown that all functional 
dependencies on the coordinates yi remain the same. 
Only expressions of coefficients in terms of deriva- 
tives of potential have changed. 

FIRST APPROXIMATION NEAR CUSP 

We assume that the origin of coordinates in eq. ([2]) 
is a cusp: 61 = 0. In this case, parameter of 
proximity is introduced by the following relations: 

yi = t^yi, 2/2 = t^m, xi = t'^xi, X2 = tx2. it can 
be shown that coordinates of image Xi (with param- 
eterization proposed above) are analytical functions 
oit. To return to initial coordinates, it is enough to 
put t = 1. We can find from formulae ([2]) , accurate 
within first order terms, that the lens equations near 
cusp caustic are 



yi = 2axi — 0X2 -1- (26x1X2 



dxfj 



t, 

Z.2 



20x1X2 + CX2 + (65;i — 3dxiX2 + fx2) ■ t 



(25) 



where = 02, 6 = 62, c= 02- 

We looked for solutions in the form: xi = xio + 
txii, X2 = X20 + tx2i- The basis for solutions con- 
struction is a cubic equation for X20: 



51 = —lla\b\b2 + 401026^ -|- 30o26i(i(j— 

- 7 (o^ - C2a) ^ - 46i (36ici + 62C2 + 3/a) a, (18) 

52 = 26i [3a^ + 5026162 - 2oi6? 

-3 (02C2 + 6id) a] . (19) 



Take notice that formula for Ji was found in [10] . Fi- 
nally, for the total magnification factor of two critical 
images, we obtained: 



^ 1 Q{y2) 

2o-V|6i| 2/2 



2 ■ 

l + Py2 + Qyi-7^ 
4 2/2. 



(20) 



CX20 - a-yxX2'd - oy2 = 0, 



(26) 



where C = ca — a?. 

Equation ([26|) has one or three real roots depend- 



yh' 



,1 . r ■ ^ ^2" O^^i/l 

mg on the sign of expression Q = ^7^' 

real root when Q > and three real roots when 
Q < 0. And explicit expressions for solutions of 
eq. ()26p are given with Cardano formulae. 

For the first coordinate in zero order approxima- 
tion, we found: 



1 

2^ 



2/1 + ax 



20J 



(27) 
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We do not present intermediate formulae for the first 
order corrections in form that repeats results of [3]. 
Instead, we give final and simplified expressions at 
once, which can be checked using substitution into 
eq. (|25p . Hence, we have: 



X2\ 



BiVixIq + B2y2X2o + Cbyf 
AaCE 

AaC^E 

Here the following notations are imposed: 



(28) 
(29) 



E = ayi- 3Cx^o> (30) 

Bi = 6aabc - aH - Aaa'^d - Qa^cd + Aa^af, (31) 

B2 = a {ba'^b - lOaad + Aa"^/) , (32) 

Ai = aa {5bc^ - Wacd + 4a V) , (33) 

A2 = a%-2aa%c+a^ {6bc^ - Wacd + Aa'^f) . (34) 

For Jacobian components J = t"^ (Jq + tJi) we 
found the following expressions: 



Ji 



Jo = -2E, 

h (3Cx^o + ayi) §2 + hx2Qyi 
CE 



(35) 
(36) 



where 



Ii =aH + a (lOad - 66c) - Aa'^f, (37) 

I2 = 16a^d-8a'^bc-a'^ {Qacd - 3bc^ + 4aV) • (38) 

The magnification factor of each image in the first 
approximation is given by the expression: 

|J| t2 |Jo + tJl| t^\Jo\ V Jo J 

While finding last equality, we took into account that 
\tJi/Jo\ < 1 (for small values of parameter t). 

RESULTS AND CONCLUSIONS 

The obtained formulae ([7|)- (|19p represent expres- 
sions of the second-order corrections for image coor- 
dinates and Jacobian near fold caustic in the case of 
general eq. ([T]). Formulae (|2U|) -(l24 p describe the to- 
tal magnification of two critical images in the second- 
order approximation with respect to proximity to the 



caustic. It is important to note that the functional 
dependence on the coordinates yi and on fitting pa- 
rameters remain the same, as in the case of k (x) = 0. 
All the differences are in expressions for the fitting 
parameters; these expressions have four additional 
constants when a continuous matter is distributed 
near the line of sight. The same situation will be with 
formulae for the magnification factor of extended 
sources [1] |3] provided that we correspondingly re- 
place coefficients P,Q,k and take into account that 
o" 7^ 1. Coefficients that are discussed in the present 
paper play a role of adjustable parameters in mod- 
elling observable light curves. Specifically, taking 
into account a continuous matter does not change 
anything in previous treatment of the strong magni- 
fication event in Q2237+0305 [,1][3]. Explicit depen- 
dencies of coefficients (I2ip -(p^ on the derivatives of 
potential $(x) will be important in case of modelling 
deflector mass distribution. 

In the last section we obtained the first-order cor- 
rections for the image coordinates and the Jacobian 
near a cusp caustic (|28p -(f38]). Some inaccuracies of 
paper [4J have been corrected, and explicit expres- 
sions of the corrections are found in terms of the 
potential expansion parameters and the roots of the 
cubic equation (j26p . 
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